Controllable coupling of superconducting flux qubits 
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We have realized controllable coupling between two three-junction flux qubits by inserting an 
additional coupler loop between them, containing three Josephson junctions. Two of these are 
shared with the qubit loops, providing strong qubit-coupler interaction. The third junction gives 
the coupler a nontrivial current-flux relation; its derivative (i.e., the susceptibility) determines the 
coupling strength J, which thus is tunable in situ via the coupler's flux bias. In the qubit regime, 
J was varied from ~45 (antiferromagnetic) to ^—55 mK (ferromagnetic); in particular, J vanishes 
for an intermediate coupler bias. Measurements on a second sample illuminate the relation between 
two-qubit tunable coupling and three-qubit behavior. 

PACS numbers: 85.25.Cp, 85.25.Dq, 03.67.Lx 



The development of Josephson qubit devices has led 
to the realization of quantum gates []], 0, 0] as well as 
two- [l], 0, [H[ and four-qubit @ coupling. For the imple- 
mentation of real quantum algorithms, several hurdles 
must still be cleared, such as increasing the number of 
qubits and their coherence times. Equally important, 
however, is coupling tunability. If the coupling strength 
can be continuously tuned between two values with op- 
posite signs, it can be naturally switched off — a great ad- 
vantage when applying two-qubit gates [7]. Moreover, in 
adiabatic quantum computing, continuously tuning the 
Hamiltonian is crucial, and both ferromagnetic (FM) and 
antiferromagnetic (AF) couplings are necessary [8[- In 
a promising group of proposals, coupling capacitances 
and inductances are replaced with effective ( "quantum" ) 
capacitances [9] and inductances EH, respectively, 
which are (sign-)tunable via their bias dependence. 

We report the realization of sign-tunable coupling 1 121 
between three- Josephson-j unction (3JJ) flux qubits [131 ]. 
These have a low charge-noise sensitivity, common to all 
flux qubits. Their small area also protects them reason- 
ably well against magnetic noise, but limits the strength 
of their AF coupling via magnetic [5[ and/or kinetic [14( 
inductance. This can be overcome by using a Joseph- 
son mutual inductance [HI , which can also be "twisted" 
for FM coupling or (in theory) current-biased for lim- 
ited tunability [u|. Our design [13] combines the above 
ideas, mediating a tunable galvanic coupling through a 
"quantum Josephson mutual inductance". The coupler 
is connected to qubits a and b via shared junctions 7 
and 9, see Fig. [TJ By changing the coupler's flux bias 
= f c &o ($o is the flux quantum), the phase difference 
across junction 8 and therefore the interaction strength J 
can be tuned. The fluxes through the coupler and qubits 




FIG. 1: Circuit design of sample 1. Junctions 123 (456) form 
qubit a (b). The large qubit junctions have areas Si, 3, 4, 6 = 
150 x 700 nm 2 , while a a , b = S 2 ,s/Si = 0.65. Junctions 789 
form the coupler (red), with areas £7,9 = 150 x 2000 nm 2 
and 5*8 = 150 x 400 nm 2 . The ratio between the coupler and 
qubit loop areas is 1 : 2. The tank circuit has inductance 
L = 90 nH and capacitance C = 470 pF, yielding a resonance 
at cj t /2tt = 20.76 MHz with quality Q = 300. The mutual 
inductances, deduced from the flux periodicity, are M a ,T = 
M b , T 98, M c ,t 51, M aM 0.85, M bM 0.17, M cM 
0.19, M a , b2 = 0.23, M b , b2 = 1.2, M c , b2 = 0.25, M a , b3 = 0.42, 
M fe , b3 = 0.39, and M c , b3 = 0.56 (all in pH). 



are controlled by bias-line currents 2bi,2,3 an d the dc com- 
ponent ibT of the coil current [181 ]. 

The system can be described by the effective pseu- 
dospin Hamiltonian 0, H, EH 

H = - E + A ^ + J (^>i a M b) , (i) 



where q is the bias on qubit z, is the corresponding 
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FIG. 2: Black: the current-flux relation I(f c ) of a coupler 
with a c = Ss/Sg = 0.2 and I c — 1 fiA. Red-dashed line: 
the coupling energy J(/ c ) obtained from Eq. J3} using this 
I(fc) and the loop currents I pa ,b found independently from 
the qubit response. Circles: experimental J(f c ) obtained from 
Fig. El 



tunnelling matrix element, ai l \ are Pauli matrices 
in the flux basis, and C = 27r/ c . 

To calculate J for a coupler with symmetric Josephson 
energies £7^9 = E, E$ = a c E, consider the potential 

f/j(0) = -£[cOS0 7 + a c COs(0 c -07-0 9 ) +COS0 9 ] 

+ ^a(0a+0 7 ) + Ubifo+fa) , (2) 

where we implemented flux quantization for small- 
inductance loops, with (j) a ^ being the qubit flux biases in 
phase units, and where the qubit energies U a ,b are already 
minimized over their internal degrees of freedom 0i 6- 
However, each qubit has two minimum states, with op- 
posite values of the persistent currents —(2e/h)U' ab = 
±^pa,6- We minimize Uj(4>) with respect to 07,9, and 
expand in E~ x . This implements the classical limit of 
the general condition that the coupler should stay in 
its ground state, following the qubits adiabatically 
To leading order, the phases obey 07 = 09 = 0, with 
sin0 = a c sin(0 c — 20). Proceeding to 0(E~ 1 ) and re- 
taining terms ex I pa Ipb, one finds [19[ 



J(<f>c) 



2el? 



(3) 



in terms of the coupler current I((f) c ) — I c sin0 with I c = 
(2e/K)E, so that 



cos(0)a c cos(0 c — 20) 
' cos(0) + 2a c cos(0 c — 20) 



(4) 



The numerator in Eq. (|3]) also occurs for magnetic cou- 
pling 11]; the denominator reflects, for finite coupler- 
loop currents, the nonlinear ity of the Josephson elements 
7 and 9. Figure [2] shows / and J for a c = 0.2. If 



a c < 1, /(0c) « a c / c sin0 c , i.e., J(0 C ) « I pa Ipb x 
h 2 a c cos c /4e 2 £ [17]. Hence, J(0 C ) > (< 0) near C = 
(7r) [20], corresponding to AF (FM) coupling. However, 
Fig. [2] already is strongly non-sinusoidal, with a larger 
maximum for FM than for AF coupling (cf. Ref. [111]). 

The qubit-coupler circuit was fabricated out of alu- 
minum, and the pancake coil out of niobium 0, 0] . Be- 
sides providing an overall dc field bias, the coil is part of 
an LC tank circuit driven at resonance, well below the 
characteristic qubit frequencies. In the Impedance Mea- 
surement Technique [2l|, one records the tank's current- 
voltage phase angle which is very sensitive to its effec- 
tive inductance: tan# = — (Q/L)x\ where Q and L are 
the tank's quality factor and "bare" inductance, respec- 
tively, and x' is the contribution to the tank inductance 
due to the qubits' reactive susceptibility. In the coherent 
regime, x' has peaks at level anticrossings 22], where a 
small bias ei ~ A* will flip the flux state of qubit i. Thus, 
these peaks demarcate the qubits' stability diagram, from 
which J can be read off [16], while their widths are cx A^. 

All measurements were performed in a dilution refrig- 
erator with a base temperature of 10 mK. Results for 
sample 1 are presented in Fig. [3] around the qubits' co- 
degeneracy point [e a = 65 = in (pQ)], for fixed [HI cou- 
pler fluxes f c = -0.08, -0.30, and -0.38. In Fig.Ek, the 
FM ordered |||) and | jj) states are pushed away from 
this point, where the |||) and |||) states dominate — a 
clear signature of AF coupling. In Fig. [3)3, the a- (ver- 
tical) and 6-traces (horizontal) are independent, demon- 
strating zero coupling. Finally, Fig. [3b is opposite to 
Fig. [3^, corresponding to FM coupling. 

Quantitatively, the state of the system (pQ) at tempera- 
ture T e ff is readily calculated; the effect on the tank flux 
follows from the mutual inductances as in Fig.[IJ and tak- 
ing the JbT-derivative yields x' • Fitting this equilibrium 
response to the data, all parameters in Eq. (pQ) as well as 
T e ff can be extracted 

d EE & From the shape of 
the single-qubit traces, one finds A ay b and I pa ,b] the shifts 
Af a ^ in these traces when they cross the co-degeneracy 
point yield J = ^Af a ^oI pa ,b- That is, here we fore- 
most measure H, not the qubit state As an exam- 
ple, Fig. [4] shows a fit to Fig. [3^, yielding T e ff ~ 70 mK, 
A a = 300 mK, J pa = 75 nA, A 6 = 55 mK, I ph = 180 nA, 
and J ~ 45 mK. The agreement between theory and ex- 
periment confirms that the system is in the qubit regime. 

The J(/ c ) thus measured (Fig. [2]) agrees with Eq. (j3]) 
for the design value a c = 0.2, and I c = 1 fiA for the 
coupler. Note that J(/ c =— 0.30) = already implies 
a c ~ 0.2 by Eqs. ([3]) and I c = 1 \i A is consistent 
with J c w (Sg/S^Ipb/V 1 ~ V 4 ^ 2 ~ °- 8 expected 
from I v b (*56, £9 are the respective junction areas). 

A sample 2 was fabricated with a c = 0.5, i.e., a larger 
junction 8 and hence potentially a stronger coupling. 
Figure [5] presents results for tan#(/bT, ^bi)- In Fig. [5^, 
Ib2 — 0. The qubits are at co-degeneracy for f c = 0.237. 
A theoretical fit as for sample 1 gives I pa = 130 nA, 
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FIG. 3: — tan 0(f a , fb) for sample 1 at coupler bias f c = —0.08, —0.30, and —0.38, with a manifest change in coupling sign. A 
theory fit as in Fig. |4]yields couplings J = 45, 0, and —55 mK. The excess response in the |ff) quadrant for f c = —0.38 is due 
to the coupler. 



I P b = 70 nA, A a « A b « 70 mK, and J « 85 mK. In 
Fig. [5b, a bias current Jb2 = 450 /iA is applied. As a 
result, the co-degeneracy point is shifted to f c = 0.347. 
Compared with Fig. [5^, the coupling strength is reduced 
to J = 85 x (0.30/0.58) « 44 mK. 

For a c = 0.5, the negative-slope portion of I(f c ) (cf. 
Fig. [2]) is very narrow, so that FM coupling should only 
occur for \f c — -|| <C 1. Hence, in sample 2 only the (AF) 
coupling strength could be varied, not its sign. In fact, 
the coupler is on the very boundary of the hysteretic 
regime [13|- Hence, for f c « ~, the energy gap above the 
ground state will be very small and the adiabaticity con- 
dition mentioned above Eq. (|3|) breaks down, regardless 
of the exact value of a c — ^. Indeed, in this case we rather 
observe three-qubit behavior, with the coupler's own anti- 
crossing characterized by A c = 90 mK and 7 pc = 300 nA 
(Fig. [5t). This is fully consistent with our interpretation " 400 
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FIG. 4: Theoretical fit for Fig. [3^l. The extracted parameters 
are T eff 70 mK, A a 300 mK, I pa 75 nA, A b 55 mK, 
I pb = 180 nA, and J(-0.08) « 45 mK. 



above, since the operating regimes are different. 

In conclusion, we have for the first time demon- 
strated sign-tunable Josephson coupling between two 
three-junction flux qubits, in the quantum regime. At 
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FIG. 5: — tan^(/bT, Ihi) for sample 2. (a) ib2 = 0. At co- 
degeneracy, the coupler flux is f c = 0.237, and the induced 
coupling J w 85 mK. (b) and (c) I h2 = 450 fiA. (b) Close- 
up of the co-degeneracy point [dashed box in (c)], where 
f c = 0.347 and J « 44 mK. (c) Overview and stability di- 
agram. At bias points outside the boxed area, the coupler 
itself can become bistable, and contribute its own anticross- 
ing to the impedance response. The arrow sizes denote the 
relative magnitudes of the loop currents in qubit a, coupler, 
and qubit b respectively. 
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T e ff w 70 mK, the coupling strength J was changed from 
+45 (antiferromagnetic) to —55 mK (ferromagnetic). At 
an intermediate coupler bias, J vanishes, thereby real- 
izing the elusive superconducting switch. These results 
represent considerable progress towards solid-state quan- 
tum computing in general. The present low-frequency 
mode of operation is particularly attractive for adiabatic 
quantum computing: control of J > is necessary to op- 
erate the computer, and sufficiently strong ferromagnetic 
coupling (J < 0) allows one to create dummy qubits, as 
used in the scalable architecture of Ref. While our 
measurements are essentially equilibrium, the design of 
Fig. [D is also relevant in the ac domain, where the cou- 
pling can be controlled by a resonant rf signal [231 ]. 
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Note added. — Recently, we learned of the work of Hime 
et al [24| implementing the controllable-coupling pro- 
posal described in Ref. [101]. 
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